**********************************************************************************************
set more off
log using table5.log, replace text


/*
*instructions for outreg;
local instruct_ols "ctitle(OLS) dec(4)  auto(2) symbol(***,**,*) se"; 
local instruct_iv "ctitle(IV) dec(4)  auto(2) symbol(***,**,*) se addstat(First stage F, e(widstat))"; 
local instruct_iv_overid "ctitle(IV) dec(4)  auto(2) symbol(***,**,*)  se addstat(Overid,e(jp),First stage F, e(widstat))"; 
*/


********************************************************************; 
***     TABLE 5 pop growth 1990_2010 -- IV excluding the West          
********************************************************************; 

*stripped down specification
ivreg2 D_censuspop9010_cp (all_road_2010_rays=road_1962_rays) if west==0 & D0==1,  first ffirst robust cluster(province05)
outreg2 using "table5", excel  dec(4)  auto(2) symbol(***,**,*)  se addstat(First stage F, e(widstat))  replace
*base specification
ivreg2 D_censuspop9010_cp larea_cc90 larea_pf05 province_capitalPlus lcensuspop1982_pref frac_highedu1982_pref share_emp_man1982_pref (all_road_2010_rays D_censuspop9010_pref=road_1962_rays ruralMigp00) if west==0 & D0==1,  first ffirst robust cluster(province05)
outreg2 using "table5", excel  dec(4)  auto(2) symbol(***,**,*)  se addstat(First stage F, e(widstat))  append
*add km roads outside CC
ivreg2 D_censuspop9010_cp larea_cc90 larea_pf05 province_capitalPlus lcensuspop1982_pref frac_highedu1982_pref share_emp_man1982_pref (all_road_2010_rays D_censuspop9010_pref lall_road_2010_km_pfcp=road_1962_rays ruralMigp00 lroad_1962_km_pfcp) if west==0 & D0==1, first ffirst robust cluster(province05)
outreg2 using "table5", excel  dec(4)  auto(2) symbol(***,**,*)  se addstat(First stage F, e(widstat)) append
*add radial rr instead
ivreg2 D_censuspop9010_cp larea_cc90 larea_pf05 province_capitalPlus lcensuspop1982_pref frac_highedu1982_pref share_emp_man1982_pref (all_road_2010_rays D_censuspop9010_pref rail_2010_rays=road_1962_rays ruralMigp00 rail_1962_rays) if west==0 & D0==1, first ffirst robust cluster(province05)
outreg2 using "table5", excel  dec(4)  auto(2) symbol(***,**,*)  se addstat(First stage F, e(widstat)) append
*ring-road only
ivreg2 D_censuspop9010_cp larea_cc90 larea_pf05 province_capitalPlus lcensuspop1982_pref frac_highedu1982_pref share_emp_man1982_pref (po_s_all_road_2010_ringxda D_censuspop9010_pref = po_s_road_1962_ringxda ruralMigp00 ) if west==0 & D0==1, first ffirst robust cluster(province05)
outreg2 using "table5", excel  dec(4)  auto(2) symbol(***,**,*)  se addstat(First stage F, e(widstat)) append
*radial + ring-roads
ivreg2 D_censuspop9010_cp larea_cc90 larea_pf05 province_capitalPlus lcensuspop1982_pref frac_highedu1982_pref share_emp_man1982_pref (all_road_2010_rays po_s_all_road_2010_ringxda D_censuspop9010_pref = road_1962_rays po_s_road_1962_ringxda ruralMigp00 ) if west==0 & D0==1, first ffirst robust cluster(province05)
outreg2 using "table5", excel  dec(4)  auto(2) symbol(***,**,*)  se addstat(First stage F, e(widstat)) append



log close

